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Abstract 

This paper aims to give explicitly all non-linear vacuum solutions 
to our non-linear field equations [1], and to define in a coordinate free 
manner the important subclass of non-linear solutions, which we call 
almost photon-like. By means of a correct introduction of the local and 
integral intrinsic angular momentums of these solutions, we saparate 
the photon-like solutions through the requirement their integral in- 
trinsic angular momentums to be equal to the Planck's constant h. 
Finally, we consider such solutions, moving radially to or from a given 
center, using standard spherical coordinates. 
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1 Explicit non-linear vacuum solutions 

As it was shown in [1] with every nonlinear solution F of our nonlinear 
equations (we use all notations from [1]) 

SF A *F = 0, 5*FA**F = 0, 5 * F A *F - 5F A F = 0. (1) 

a class of F-adapted coordinate systems is associated, such that F and *F 
acquire the form respectively 

F = eudx Adz + udx A d^ + epdy Adz + pdy A d£ 

*F = —pdx Adz — epdx Ad^ + udy Adz + eudy A d^. 

Since we look for non-linear solutions of (1), we substitute these F and *F 
in (1) and after some elementary calculations we obtain 

SF = (v,£ - eu z )dx + (p£ - ep z )dy + e(u x + p y )dz + (u x + p y )d£, 

5* F = -e(p£ - ep z )dx + e(u^ - ep z )dy - (p x - u y )dz - (p x - u y )d£, 
F^{5FYdx» = (*F)^(5 * F) v dx v = 
= £ \p(pn ~ £Vz) + u(it{: - eu z )\ dz + \p(p£ - ep z ) + u(u£ - eu z )\ d£, 
(SF) 2 = (5 * Ff = -( Ui - eu z ) 2 - - ep z ) 2 
A simple direct calculation shows, that the equation 

5*FA*F-5FAF = 

is identically fulfilled for any such F and *F with arbitrary u and p. We 
obtain that our equations reduce to only 1 equation, namely 

p(Pt ~ £Pz) + u(u^ - eu z ) = i [(u 2 + p 2 )(: - e(u 2 + p 2 ) z ] = 0. (2) 

The obvious solution to this equation is 

u 2 +p 2 = (P 2 (x,y,£ + ez). (3) 

The solution obtained shows that the equations impose some limitations 
only on the amplitude function and the phase function f is arbitrary exept 
that it is bounded: \ip\ < 1. The amplitude is a running wave along the 
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specially chosen coordinate z, which is common for all F-adapted coordinate 
systems. Considered as a function of the spatial coordinates, the amplitude 
is arbitrary, so it can be chosen spatially finite. The time-evolution does not 
affect the initial form of 0, so it will stay the same in time. This shows, that 
among the nonlinear solutions of our equations there are (3+1) soliton-like 
solutions. The spatial structure is determined by the initial condition, while 
the phase function <p can be used to define internal dynamics of the solution. 
Recalling from [1] the substitutions 

U = (j).ip, p = <j)yj\ - V? 2 , 

and the equality \A\ — 0, we get 



For the induced pseudoorthonormal bases (1-forms and vector fields) we find 
A = pdx + \J\ — p 2 dy, eA* = — tp 2 dx + ipdy, R = — dz, S = d£, 



Hence, the nonlinear solutions in canonical coordinates are parametrized 
by one function of 3 parameters and one bounded function of 4 parameters. 
Therefore, the separation of various subclasses of nonlinear solutions is made 
by imposing additional conditions on these two functions. Further we are 
going to separate a subclass of solutions , the integral properties of which 
reflect well enough the well known from the experiment integral properties 
and characteristics of the free photons. These solutions will be called photon- 
like and will be separated through imposing additional requirements on and 
ip in a coordinate-free manner. 



2 Almost photon-like solutions 

We note first, that we have three invariant quantities at hand: 0, ip and 
L. The amplitude function is to be determined by the initial conditions, 
which have to be finite. So, we may impose additional conditions on L and p. 
These conditions have to express some intra- consistency among the various 
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characteristics of the solution. The idea, what kind of intra-consistency to 
use, comes from the observation that the amplitude function (ft is a first 
integral of the vector field V, i.e. 

~ £ ~d~z + di) ^ = ~ e ~d~ z ^ x ' + £z) + d~^ x ' + ez) = °' 

We want to extend this available consistency between V and 0, so we shall 
require the two functions p and L to be first integrals of some of the available 
F-generated vector fields. Explicitly, we require the following: 

1°. The phase function ip is a first integral of the three vector fields A, A* 
and R: A(y?) = 0, A*((p) = 0, R(<p) = 0. 

2°. The scale factor L is a non-zero finite first integral of the vector field 
S: S(L) = 0. 

The requirement R(y?) = just means that in these coordinates p> does 
not depend on the coordinate z. The two other equations of 1° define the 
following system of differential equations for p>: 



dtp r ~ 2 dp> r -dp dip 

Noticing that the matrix 



-ip -Vi - tp 2 

has non-zero determinant, we conclude that the only solution of the above 
system is the zero-solution: 

dp = dp =Q 
dx dy 

We obtain that in the coordinates used the phase function p depends only 
on £. Therefore, in view of (4), for the scale factor L we get 

m\ 
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Now, the requirement 2°, which in these coordinates reads 

s(£) « » vr^ = „, 

just means that the scale factor L is a pure constant: L = const. In this way 
we obtain the differential equation 

The obvious solution to this equation reads 

= cos \ k t + const J , (6) 



where k — ±1. It worths to note that the naturally arising characteristic 
frequency according to the equation 

" = !■ < 7) 

has nothing to do with the concept of frequency in CED. In fact, the quantity 
L can not be defined in Maxwell's theory. 

Finally (recalling [1]) we note, that the so obtained phase function </?(£) 
leads to the following. The 2-form tr(TZ°), where TZ° is the matrix of 2-forms, 
formed similarly to the matrix 1Z, but using the basis (A, eA*, R, S) instead 
of the basis (A, eA*, R, S), is closed. In fact, 

tr(R°) = (pdx Ad£ + (pdy A d£ - dy A dz + dz A d£ 

and since f = <£>(£), we get dtr{lZ Q ) = 0. Note also that the above explicit 
form of tr(TZ°) allows to define the phase function by 



up = 



\tr(R { 



0\\2 



2 

This class of solutions we call almost photon-like. 

Remark. If one of the two functions u and p, for example p, is equal to zero: 
p — 0, then, formally, we again have a solution, which may be called linearly 
polarized by obvious reasons. Clearly, the phase function of such solutions 
will be constant: <p = const, so, the corresponding scale factor becomes 
infinitely large: L — > oo, therefore, condition 2° is not satisfied. The reason 
for this is, that at p = the function u becomes a running wave and we get 
\SF\ = \5 * F\ =0, so the scale factor can not be defined by the relation 
L = \A\/\SF\. 
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3 Intrinsic angular momentum (helicity) 
and photon-like solutions 

The problem for describing the intrinsic angular momentum (I AM), or in 
short helicity, spin of the photon is of fundamental importance in modern 
physics, therefore, we shall pay a special attention to it. In particular, we 
are going to consider two approaches for its mathematical description. But 
first, some preceding comments. 

First of all, there is no any doubt that every free photon carries such an 
intrinsic angular momentum. Since the angular momentum is a conserved 
quantity, the existence of the photon's intrinsic angular momentum can be 
easily established and, in fact, its presence has been experimentally proved 
by an immediate observation of its mechanical action and its value has been 
numerically measured. Assuming this is so, we have to understand its origin, 
nature and its entire meaning for the existence and outer relations of those 
natural entities, called shortly photons somewhere in the first quarter of this 
century. 

So, we begin with the assumption: every free photon carries an intrin- 
sic angular momentum with integral value equal to the Planck's constant h. 
According to our understanding, the photon's IAM comes from an intrinsic 
periodic process. This point of view undoubtedly leads to the notion, that 
photons are not point-like structureless objects, they have a structure, i.e. 
they are extended objects. In fact, according to one of the basic principles 
of physics all free objects move as a whole uniformly. So, if the photon is 
a point-like object any characteristic of a periodic process, e.g. frequency, 
should come from an outside force field, i.e. it can not be free: a free point- 
like (structureless) object can not have the characteristic frequency. 

This simple, but true, conclusion sets the theoretical physics of the first 
quarter of this century faced with a serious dilemma: to keep the notion of 
structurelessness and to associate in a formal way the characteristic frequency 
to the microobjects, or to leave off the notion of structuelesness, to assume 
the notion of extendedness and availability of intrinsically occurring periodic 
process and to build corresponding integral characteristics, determined by 
this periodic process. A look back in time shows that the majority of those 
days physicists had adopted the first approach, which has brought up to life 
quantum mechanics as a computing method , and the dualistic-probabilistic 
interpretation as a philosophical conception. If we set aside the widespread 
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and intrinsically controversial idea that all microobjects are at the same time 
(point-like) particles and (infinite) waves, and look impartially, in a fair- 
minded way, at the quantum mechanical wave function for a free particle, we 
see that the only positive consequence of its introduction is the legalization of 
frequency, as an inherent characteristic of the microobject. In fact, the prob- 
abilistic interpretation of the quantum mechanical wave function for a free 
object, obtained as a solution of the free Schroedinger equation, is impossible 
since its square is not an integrable quantity (the integral is infinite). The 
frequency is really needed not because of the dualistic-probabilistic nature 
of microobjects, it is needed because the Planck's relation E = hu turns 
out to be universally true in microphysics, so there is no way to avoid the 
introduction of frequency. The question is, if the introduction of frequency 
necessarily requires some (linear) "wave equation" and the simple complex 
exponentials of the kind const. exp[i(k.r — id)], i.e. running waves, as "free 
solutions". Our answer to this question is "no". The classical wave is some- 
thing much richer and a much more engaging concept, so it hardly worths to 
use it just because of the attribute of frequency. In our opinion, it suffices to 
have a periodic process at hand. 

These considerations made us turn to the soliton-like objects, they present 
the two features of the microobjects (localized spatial extendedness and time- 
periodicity), simultaneously, and, therefore, seem to be more adequate theo- 
retical models for those microobjects, obeying the Planck's relation E = hu. 
Of course, if we are interested only in the behaviour of the microobject as a 
whole, we can use the point-like notion, but any attempt to give a meaning 
of its integral characteristics without looking for their origin in the consis- 
tent intrinsic dynamics and structure, in our opinion, is not a perspective 
theoretical idea. And the "stumbling point" of such an approach is just the 
availability of an intrinsic mechanical angular momentum, which can not be 
understood as an attribute of a free structureless object. 

Having in view the above considerations, we are going to consider two 
ways to introduce and define the intrinsic angular momentum as a local 
quantity and to obtain, by integration, its integral value. So, these two ap- 
proaches will be of use only for the spatially finite nonlinear solutions of 
our equations. The both approaches introduce in different ways 3-tensors 
(2-covariant and 1-contravariant). Although these two 3-tensors are built of 
quantities, connected in a definite way with the field F, their nature is quite 
different. The first approach is based on an appropriate tensor generaliza- 
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tion of the classical Poynting vector. The second approach makes use of the 
concept of torsion, connected with the field F, considered as 1-covariant and 
1-contravariant tensor. The first approach is pure algebraic, while the sec- 
ond one uses derivatives of F^. The spatially finite nature of the solutions 
F allows to build corresponding integral conserved quantities, naturally in- 
terpreted as angular momentum. The scale factor L appears as a multiple, 
so these quantities go to infinity for all linear (i.e. for Maxwell's) solutions. 

In the first approach we make use of the scale factor L, the isotropic 
vector field V and the two 1-forms A and A*. By these four quantities we 
build the following 3-tensor H: 

H = k^V ®(AAA*). (8) 

The connection with the classical vector of Poynting comes through the ex- 
teriour product of A and A*, the 3-dimensional sense of which is just the 
Pointing's vector. In components we have 

HZ a = K^V»(A u A* a -A a A:). 
In our system of coordinates we get 

E = K~ c + <8> (e<f) 2 dx A dy), 

so, the only non-zero components are 

Hf 2 = -Hi, = -K-cf, Hf 2 = -Hi = ke-<P 2 . 

c c 

It is easily seen, that the divergence V^H^ a — > V M iff 2 is equal to 0. In fact, 



V,#f 2 = ^-H\ 2 + ^H' l2 = k- c [-(<j> 2 ) z + (e^ 







because <p 2 is a running wave along the coordinate z. Since the tangent 
bundle is trivial we may construct the antisymmetric 2-tensor 



/ H Ava dxdydz, 
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the constant components of which are conserved quantities. 



H 12 — — H21 — 



H^dxdyd, 



\ Z = KE-W = KEWT = KE-. 

C V 



The non-zero eigen values of H„ CT are pure imaginary and are equal to ±iWT. 
This tensor has unique non-zero invariant P(F), 



The quantity P(F) will be called Planck's invariant for the finite nonlinear 
solution F. All finite nonlinear solutions Fi, F 2 , , satisfying the condition 



where h is the Planck's constant, will be called further photon-like. The 
tensor field H will be called intrinsic angular momentum tensor and the 
tensor H will be called spin tensor or helicity tensor. The Planck's invariant 
P(F) = WT, having the physical dimension of action, will be called integral 
angular momentum, or just spin or helicity. 

The reasons to use this terminology are quite clear: the time evolution of 
the two mutually orthogonal vector fields A and A* is a rotational-advancing 
motion around and along the z-coordinate (admissible are the right and the 
left rotations: k — ±1) with the advancing velocity of c and the frequency 
of circulation v = c/L. We see the basic role of the two features of the solu- 
tions: their soliton-like character, giving finite value of all integral quantities, 
and their nonlinear character, allowing to define the scale factor L correctly. 
From this point of view the intrinsic angular momentum h of a free photon 
is far from being incomprihensible quantity, connected with the even more 
incomprihensible duality "wave-particle", and it looks as a quite normal in- 
tegral characteristic of a solution, presenting a model of our knowledge of 
the free photon. 

We proceed to the second approach to introduce I AM by recalling the 
definition of torsion of two (1,1) tensors. If G and K are 2 such tensors 





P(F 1 ) = P(F 2 ) 



h, 



G = G v „dx» ® 



d 



K = K v n dx^ ® 



d 



dx v ' 
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their torsion is defined as a 3-tensor S° = 



—S° by the equation 



S(G, K)(X, Y) = [GX, KY] + [KX, GY] + GK[X, Y] + KG[X, Y\ — 

-G[X, KY] - G[KX, Y] - K[X, GY] - K[GX, Y], 
where [, ] is the Lie-bracket of vector fields, 

d „,._., , „ d 



GX 



dx v 



GK 



dx u 



and X, Y are 2 arbitrary vector fields. If G = K, in general S(G, G) ^ and 

S(G, G)(X, Y) = 2 { [GX, GY] + GG[X, Y] - G[X, GY] - G[GX, Y]} . 

This last expression defines at every point x G M the torsion S(G, G) = Sg of 
G with respect to the 2-dimensional plain, defined by the two vectors X(x) 
and Y(x). Now we are going to compute the torsion Sp of the nonlinear 
solution F with respect to the intrinsically defined by the two unit vectors 
A and eA* 2-plain. In components we have 



(Sf) 



/IV 



pa 



dF° 



M dx a 



dF° 

pa A* 



" dx a a dx» 



8F a dF« 



a dx v 



In our coordinate system 
d 



d_ 

dy 



eA* 



d 



d 



^ dx ^ dy' 



so, 



(S F )%A^eA* u = (Sp^AhA* 2 - A 2 eA* 1 ). 
For (Sp) 12 we get 



( s f){ 2 = (Sf)1 2 = 0, (S F )l 2 = - £ ( S F)t 2 = 2e{p(uf - eu z ) - u(p^ - ep z )}. 



Remark. In our case (Sf) 12 — (S*f)uj so further we shall work with Sf 
only. 

It is easily seen that the following relation holds: A^A* 2 — A^A* 1 = 1. 
Now, for the almost photon-like solutions 

u = <p(x, y, £ + ez) cos (^ K ~j^ + const^j , p = (f)(x, y, £ + ez) sin (^ K J^ + cons t^j 
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we obtain 

(S F ) 3 12 = -e{S F )\ 2 = -2e^\ 



{S F )%A»eA* v = 



0, 0, -2e-<p 2 , 2-<, 



Since <fi 2 is a running wave along the z-coordinate, the vector field S F (A, eA*) 
has zero divergence: V„ [S F (A, eA*)] u = 0. Now we define the helicity vector 
for the solution F by 

Z F = ^-S F (A,eA*). 
2c 

Since L = const, then T, F has also zero divergence, so the integral quantity 



J (T lF ) A dxdydz 



does not depend on time and is equal to kWT. The photon-like solutions are 
separated in the same way by the condition WT = h. Here are three more 
integral expressions for the quantity WT. We form the 4-form 

1 K 
— — S A *Y, F = —(p 2 (x! 
L c 

and integrate it over the 4-volume 1Z 3 x L, the result is kWT . Besides, we 
verify easily the relations 



- / \A A A*\uj = — [ \5F A 5 * F\cu = WT. 

C JK A xL C JRA-kL 



Since we separate the photon-like solutions by the relation WT = h, the 
last expressions suggest the following interpretation of the Planck's constant 
h. Since \A A A*\ is proportional to the area of the square, defined by the 
two mutually orthogonal vectors A and eA*, the above integral sums up 
all these areas over the whole 4-volume, occupied by the solution F during 
the intrinsically determined time period T, in which the couple (A,eA*) 
completes a full rotation. The same can be said for the couple (SF, 5 * F) 
with some different factor in front of the integral. This shows quite clearly 
the "helical" origin of the full energy W = hv of the single photon. 
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4 Solutions in spherical cordinates 



The so far obtained soliton-like solutions describe objects, "coming from 
infinity" and "going to infinity". Of interest are also soliton like solutions 
"radiated" from, or "absorbed", by some central "source" and propagating 
radially from or to the center of this source. We are going to show, that our 
equations [1] admit such solutions too. We assume this central source to be 
a small ball RP with radius r , and put the origin of the coordinate system at 
the center of the source-ball. The standard spherical coordinates (r, 9, ip, £) 
will be used and all considerations will be carried out in the region out of 
the ball BP. In these coordinates we have 

ds 2 = -dr 2 - r 2 d6 2 - r 2 sin 2 9dip 2 + d£ 2 , ^\ = r 2 sin9. 

The *-operator acts in these coordinates as follows: 

*dr = r 2 sin9d9 A dip A d£ *(dr A dd A dip) = (r 2 sinO)' 1 dt; 

*d9 = —sinOdr A dip A d£ *(dr A d6 A d£) = sinOdip 

*d<p = (sin6y l dr A d6d^ *(dr A dip A dg) = -(sin6y l d6 

*d£ = r 2 sin9dr A d9dip *(d6 A dip A d£) = (r 2 sin9)~ 1 dr 

*(dr A d9) = —sinQdip A d£ *(d9 A dp) = -(r 2 sin0)~ 1 dr A d£ 
*(dr A dp) = (sin6)~ 1 d6 A d£ *(d9 A d£) = -sinOdr A dip 
*(dr A d£) = r 2 sin6d6 A dip *(dip A d£ = (sin6)~ 1 dr A d6. 

We look for solutions of the following kind: 

F = eudr A d6 + ud9 A d£ + epdr A dp + pdip A d£, (10) 

where u and p are spatially finite functions. We get 

p p 

*F = -dr A d6 + e -d0 A dS, — usinOdr A dip — esinOdip A d£. 

smO sinO 

The following relations hold: 

F A F = 2e(up - up)dr A dd A dp A d£ = 0, 

F A *F = \-u 2 sinQ + u 2 sin9 - + ) dr A d9 A dip A d£ = 0, 
\ sint) smv J 

i.e. the two invariants are equal to zero: (*F) flu F^ u = 0, F^ V F^ U = 0. 
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After some elementary computation we obtain 

5FAF = 5*FA*F = e[u (ep r + p^) -p (eu r . + u^)\ dr Ad6 A dip+ 

+ [u (eu r + u^) — u (ep r + p^)] d6 A dip A d^ 

p (ep r + pz) 



F A *dF = e 



u (eu r + u^) sinO + 



— e 



u (eu r + u^) sin9 + 



sinO 
p (ep r + p 6 ) 



(*F) A *d * F = 



u (eu r + u^) sinO + 



sinO 

p (ep r + pz) 



sin9 



dr A d9 A d<p- 
dO A dcf) A , 

dr Ad9 A dip- 



u {eu r + u^) sm# + 



sinO 



dO A d<j> A d£. 



So, the two functions u and p have to satisfy the equation 

t , \ • a i P( £ Pr+Pt) n 
w + w^J sznu H — — = 0, 



sinO 



which is equivalent to the equation 

,2 \ 



u sinO 



sinO 



+ e \u sinO + 



P 



sind 



(11) 



(12) 



The general solution of this equation is 



u 2 sin9 + 



p 



sinO 



<t> 2 {£-er, 



(13) 



For the non-zero components of the energy-momentum tensor we obtain 



Qi ~ Qi ~ Qi ~ Qi 



l 



4irr 2 sin6 



u 2 sin9 + 



P 



sinO 



(14) 



It is seen that the energy density is not exactly a running wave but when we 
integrate to get the integral energy, the integrand is exactly a running wave: 

,2 



W 



— I * (QtA) = — I (u 2 sin6 + -^—Adr Ad9 A d<f>. 



Since the functions u and p are spatially finite, the integral energy W is finite, 
and from the explicit form of the energy-momentum tensor it follows the well 
known relation between the integral energy and momentum: W 2 — c 2 p 2 = 0. 



13 



References: 

1. Donev,S., Tashkova,M., Extended Electrodynamics: II. Properties 
and Invariant Characteristics of the Non-linear Vacuum Solutions, submitted 
for publication. 



14 



